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Abstract 

Finding a common factor of two multivariate polynomials with approximate coefficients is 
a problem in symbolic-numeric computing. Taking a tropical view on this problem leads to 
efficient preprocessing techniques, applying polyhedral methods on the exact exponents with 
numerical techniques on the approximate coefficients. With Maple we will illustrate our use 
of tropical algebraic geometry. 
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1 Introduction 

> 

00 . 

Tropical algebraic geometry is a relatively new language to study skeletons of algebraic varieties. 
Introductions to tropical algebraic geometry are in [32] and [48 1 Chapter 9]. Computational 
aspects are addressed in [6] and [51] . One goal of this paper is to explain some new words of this 
language, and to show how a general purpose computer algebra system like Maple is useful to 
explore and illustrate tropical algebraic geometry. For software dedicated to tropical geometry, 
we refer to Gfan [221 [23], a SINGULAR library [23], and Trim [39]. 

The roots of tropical algebraic geometry run as deep as the work of Puiseux [40J and 0s- 
^ trowski [36], therefore our focus is on answering a practical question in computer algebra: when 

do two polynomials have a common factor? Viewing this question in tropical algebraic geometry 
leads to a symbolic-numeric algorithm. In particular, we will say that tropisms give the germs 
to grow the tentacles of the common amoeba. The paper is structured in four parts, each part 
explaining one of the key concepts of the tropical sentence. 

Our perspective on tropical algebraic geometry originates from polyhedral homotopies [20J, 
|31j . |54j to solve polynomial systems implementing Bernshtem's first theorem [3J. Another related 
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approach that led to tropical mathematics is idempotent analysis. In [UJ, a Maple package is 
presented for a tropical calculus with application to differential boundary value problems. 

Related work on our problem concerns the factorization of sparse polynomials via Newton 
polytopes PQ, |17j : approximate factorization [II] . [TO] . [16J, [28], [35], and the GCD of polynomials 
with approximate coefficients [59]. The polynomial absolute factorization is also addressed in [8] 
and the lectures in [9] offer a very good overview. Criteria based on polytopes for the irreducibility 
of polynomials date back Ostrowski [36]. In this paper we restrict our examples to polynomials 
in two variables and refer to polygons instead of polytopes. The terminology extends to general 
dimensions and polytopes, see [SO] . 

That two polynomials with approximate coefficients have a common factor is quite an ex- 
ceptional situation. Therefore it is important to have efficient preprocessing criteria to decide 
quickly. The preprocessing method we develop in this paper attempts to build a Puiseux expan- 
sion starting at a common root at infinity. To determine whether a root at infinity is isolated 
or not we apply the Newton-Puiseux method, extending the proof outlined by Robert Walker 
in [56] , see also |14j . towards Joseph Maurer's general method [33] for space curves. A more 
algorithmic method than [33] is given in [2] along with an implementation in CoCoA. General 
fractional power series solutions are described in [34] . See [25], [26] and [39] for recent symbolic al- 
gorithms, and [37], [38] for a symbolic-numeric approach. The complexity for computing Puiseux 
expansions for plane curves is polynomial ]57] in the degrees. As an alternative to Puiseux series, 
extended Hensel series have been developed, see |44] for a survey. Good numerical convergence 
is reported in [2~T] . 

We show that via suitable coordinate transformations, the problem of deciding whether there is 
a common factor is reduced to univariate root finding, with the univariate polynomials supported 
on edges of the Newton polygons of the given equations. Also in the computation of the second 
term of the Puiseux series expansion, we do not need to utilize all coefficients of the given 
polynomials. In the worst case, the cost of deciding whether is a factor is a cubic polynomial in 
the number of monomials of the given polynomials. 

Certificates for the existence of a common factor consist of exact and approximate data: the 
exponents and coefficients of the first two terms of a Puiseux series expansion of the factor at 
a common root at infinity. The leading exponents of the Puiseux expansion form a so-called 
tropism [33] . The coefficients are numerical solutions of overdetermined systems. If a more 
explicit form of the common factor is required, more terms in the Puiseux expansion can be 
computed up to precision needed for the application of sparse interpolation techniques, see [13] . 
[U, [27], and [30]. 

The ConvexHull and subs commands of Maple are very valuable in implementing an interac- 
tive prototype of the preprocessing algorithm. For explaining the intuition behind the algorithm, 
we start illustrating amoebas and the tentacles. Once we provide an abstraction for the tentacles 
we give an outline of the algorithm and sketch its cost. 

Acknowledgements. This paper is based on the talk the second author gave at MICA 2008 
- Milestones in Computer Algebra, a conference in honour of Keith Geddes' 60th Birthday, 
Stonehaven Bay, Trinidad and Tobago, 1-3 May. We thank the organizers for this wonderful and 
inspiring conference for their invitation. 
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2 Amoebas 



Looking at the asymptotics of varieties gives a natural explanation for the Newton polygon. This 
polygon will provide a first classification of the approximate coefficients of the given polynomials. 
This means that at first we may ignore coefficients of monomials whose exponents lie in the 
interior of the Newton polygon. 



2.1 asymptotics of varieties 

Our input data are polynomials in two variables x and y. The set of values for x and y that make 
the polynomials zero is called a variety. Varieties are the main objects in algebraic geometry. In 
1971, G.M. Bergman [3] considered logarithms of varieties. In tropical algebraic geometry, we 
look at the asymptotics of varieties. 

log : C* x C* -> RxR 

(x,y) i-» (log(|x|),log(|y|)) " 

Because the logarithm is undefined at zero, we exclude the origin restricting the domain of our 
polynomials to the torus (C*) 2 , C* = C\{0}. Following [18], we arrive at our first new word [55] . 

Definition 2.1 (Gel'fand, Kapranov, and Zelevinsky 1994) The amoeba of a variety is its 
image under the log map. 



To see what amoebas look like, we use the plotting capabilities of Maple. We use polar 
coordinates (see |50| for more on plotting amoebas) to express a linear variety: 



t := —x H — y 
J 2 5 



1 = A :-- 
In Figured] we see the result of a Maple plot. 
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> f := l/2*x + l/5*y - 1: 

> s : = solve (f ,y) : 

> L := map(log,map(abs, [x,s] )) : 

> A := subs(x=r*exp(I*theta) ,L) ; 

> Ap := seq(plot([op(subs(theta=k*Pi/200,A)) , 

r=-100. . 100] ,thickness=6) ,k=0. .99) : 

> plots [display] (Ap,axes=none) ; 



Figure 1: The amoeba of a linear polynomial, with all Maple commands at the right. 
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2.2 compactifying amoebas lead to Newton poly topes 

We compactify the amoeba of / _1 (0) by taking lines perpendicular to the tentacles. As each line 
cuts the plane in half, we keep those halves of the plane where the amoeba lives. The intersection 
of all half planes defines a polygon. 

The resulting polygon is the Newton polygon of /. For the amoeba in Figure its compact- 
ification is shown in Figure [2j On Figure [2] we recognize the shape of the triangle, the Newton 
polygon of a linear polynomial. In general, there is a map [46] that sends every point in the 
variety to the interior of the Newton polygon of the defining polynomial equation. 

12 - 
10- 

8 - 

6- 



-2 J 2 4 6 (3 10 



Figure 2: The compactification of the amoeba: the edges of the Newton polygon (displayed at 
the right) are perpendicular to the tentacles of the amoeba. 

This geometric derivation of the Newton polygon coincides with the more formal definition. 

Definition 2.2 For f(x,y) = Cijx l y^ , Cjj £ C*. A is the support of /. The convex hull of 

A is the Newton polygon. 

The Newton polygon models the sparse structure of a polynomial. Most polynomials arising 
in practical applications have few monomials with nonzero coefficients and are called sparse. 
The Newton polygon assigns additional significance to the coefficients. Coefficients associated 
to monomials whose exponents span a vertex of the Newton polygon are more important than 
coefficients whose exponents lie in the interior of the Newton polygon. 

Plotting amoebas is actually computationally quite involved - the use of homotopy continu- 
ation methods [52] is suggested in [50]. See [35] for more about amoebas. We will see that the 
asymptotics of the amoebas will lead to a natural reduction of our problem to smaller polynomials 
in one variable. 

3 Tentacles 

The tentacles of the amoeba stretch out to infinity and are represented by the inner normals, 
perpendicular to the edges of the Newton polygon. 



(0,1) 




(0,0) (1,0) 
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3.1 the direction of the tentacles towards infinity 

Our problem may be stated as follows: Given two polynomials in two variables with approximate 
complex coefficients, is there a common factor? 

Looking at the problem from a tropical point of view, we first have the amoeba of the common 
factor in mind. The tentacles of that amoeba stretch out to infinity and these tentacles are 
represented by inner normals, perpendicular to the lines at infinity corresponding to the edges of 
the Newton polygon. Because the factor is common to both polynomials, the normals must be 
common to both polygons. So if there is a factor, there must be at least one pair of edges with 
the same inner normal vector. Such inner normal vector is a tropism, defined below. 

The tropical view will lead to solving the problem first at infinity, providing an efficient 
preprocessing criterion. We first formalize the geometric idea in a proposition. 

Proposition 3.1 Let f and g be two polynomials. If the amoebas of f and g have no tentacle 
stretching out to infinity in the same direction, then f and g have no common factor. 

Verifying the conditions of Proposition 13.11 seems nontrivial at first. 

3.2 normal fans and tropicalization 

The inward pointing normal vectors to the edges represent the tentacles of the amoeba. Consider 
for example 

/ : = x 3 y + x 2 y 3 + x 5 y 3 + x 4 y 5 + x 2 y 7 + x 3 y 7 . (3) 
In Figure [3] we show the Newton polygon of / and its normal fan. 




o A — ■ — i — ■ — i — ' — i — ' — i — ■ — i ▼ 

12 3 4 5 

Figure 3: The Newton polygon and its normal fan. 

The collection of inner normals to the edges of the Newton polygon forms a tropicalization 
of /, denoted by Trop(/). To formalize this notion, we introduce the following definitions. 
Exponents and direction vectors are related through duality via the inner product. 
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Definition 3.2 The inner product is 

<•>•) 



Z 2 xZ 2 
((i,j),(u,v)) 



z 

fit + jf . 



(4) 



Given a vector (u, u), (•, (u,v)) ranks the points For (u,v) = (1, 1), we have the usual 

degree of x l y 3 . So the direction of the tentacles are grading the points in the support. For 
example, in Figured we look at the support in the direction (—1, +!)• 




-1 x i + 


(+1) x j 


= +5 


-1 x i + 


(+1) x j 
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(+1) x j 


= -1 


-1 x i + 


(+1) x j 


= -2 



Figure 4: Grading the points in the support along (— 1, +1). 



The degree of x l y 3 in the direction (u,v) is the value of the inner product (u,v)). In 

Maple we compute weighted degrees as follows: 

Groebner [WeightedDegree] (f , [-1 , +1] , [x , y] ) ; 

The connection between Grobner bases and Newton polytopes is the subject of [IT]- This 
grading leads to homogeneous coordinates, see [I~2] . 

We arrive at a tropicalization of a polynomial via the normal fan to the Newton polygon of 
the polynomial. 

Definition 3.3 Let P be the Newton polygon of /. The inner product is denoted by (•, •). The 

normal cone to a vertex p of P is 

{ vel 2 \{0} | (p,v) =min(q,v) }. (5) 

The normal cone to an edge spanned by pi and p2 is 

{ v G R 2 \ {0} | < Plj v) = (p 2 ,v) = min(q, v) }. (6) 

qeP 

The normal fan of P is the collection of all normal cones to vertices and edges of P. All normal 
cones to the edges of P define a tropicalization of f, denoted by Trop(/). 

We speak of a tropicalization (a instead of the) because different regular triangulations of the 
Newton polygon will lead to different normal fans. We can identify the tropicalization by the 
secondary fan [18] but for our problem we do not need it. 
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4 Tropisms 



The tropical view will lead to an efficient preprocessing stage to determine whether two polyno- 
mials have a common factor. 

4.1 turning the varieties in a particular direction 

The answer to our original question Do two polynomials have a common factor? first depends 
on the relative position of the Newton polygons. We compute tropicalizations of the polynomials 
and obtain an efficient preprocessing step independent of the coefficients. 

We first want to exclude the situations where there is no common factor, already implied by 
the Newton polygons in relative general position. This is a direct consequence of Bernshtem's 
second theorem [3]. For completeness, we state this theorem here for Newton polygons. 

Theorem 4.1 Let f and g be two polynomials in x and y. If Trop(/) n Trop(g) = then the 
system f(x,y) = = g(x,y) has no solutions at infinity. 

Now we can make Proposition 13.11 effective: 

Proposition 4.2 If for two polynomials f and g: Trop(/) n Trop(g) = 0, then f and g have no 
common factor. 

Proof. By Theorem 14.11 Trop(/) n Trop(g) = implies there is no common root at infinity. But 
if / and g would have a common factor, they would have a common root at infinity as well. □ 

Consider for example the tropicalizations of two random polynomials of degree 15, shown in 
Figure 




Figure 5: The first two pictures from the left represent the normal fans of two polynomials. By 
superposition of the fans at the far right we see there are no common directions. Therefore, for 
all nonzero coefficients, the polynomials can have no common factor. 
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In our second example we generated a factor of degree 5 and multiplied with two random 
polynomials of degree 10. Tropicalization of the factor and the two polynomials are shown in 
Figure [H 





/ 



Figure 6: The normal fan at the left is the normal fan of the factor common to two polynomials / 
and g whose normal fans are displayed in the middle and at the right. We recognize the fan at 
the left as a part of the other fans. 

A dictionary definition of a tropism is the turning of all or part of an organism in a particular 
direction in response to an external stimulus. Tropisms were introduced mathematically in 1980 
by Joseph Maurer [33] who generalized Puiseux expansions for space curves. We adapt his 
definition for use to our problem. 

Definition 4.3 Let P and Q be Newton polygons of / and g. A tropism is a vector perpendicular 
to one edge of P and one edge of Q. 

Using the general terminology of [60J, tropisms correspond to the one dimensional cones in 
the common refinement of the normal fans of the polygons. 

4.2 certificates for numerical computations 

Tropisms are important because they give a first exact certificate for the existence of a common 
factor. Selecting those monomials which span the edges picked out by the tropism defines a 
polynomial system which admit a solution in (C*) 2 . 

Definition 4.4 Let (u, v) be a direction vector. Consider / = N CijX l y 3 . The initial form of 



/ in the direction (u, v) is 






(7) 



{(i,j),(u,v)) ■■ 



= m 



where m = min{ (u, v)) \ G A }. 
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The direction (u, v) is the normal vector to the line ui + vj = m which contains the edge of the 
Newton polygon of /. This edge is the Newton polygon of in^ u ^(f). 

The terminology of initial forms corresponds to the Grobner basics [47j . In [43 j . mr uv \(f) is 
called an initial term polynomial. We call a tuple of initial forms an initial form system. Initial 
form systems correspond to truncated systems in [7] and |29j. 

The factor common to the / and g generated above is 

r := 2xy + x 2 y + 9xy 2 + 7x 3 y + x 4 y + 9x 3 y 2 , (8) 
In Figure[7]we show the initial forms associated with the two polynomials defined by the tropisms. 

Investigate 4 directions, take (1,0): 
m (i,o)W = 2xy + 9xy 2 

Initial forms of / and g: 
m (i,o)(/) = 55xy 6 + 10xy 5 + 45xy 7 
' m (i,o){g) = 10xy 6 + 45xy 7 

Figure 7: The normal fan of the common factor and the initial form systems corresponding to 
the direction (1,0). 

Because the tropism is a standard basis vector (1,0), the initial form system it determines 
consists of two polynomials in one variable: 

f in (li0) (/) = x{5y 5 (y + l)(2 + 9y)) =0 

\ m {lfi) (g) = x{5y 5 (2 + 9y)) =0 { > 

and then y = —2/9 represents the common root at infinity, In general the common root at infinity 
will be an approximate root and with a-theory [5] we can bound the radius of convergence for 
Newton's method. In addition to the first certificate, the exact tropism, the common root at 
infinity is the second approximate certificate. 

For general tropisms, not equal to basis vectors, we perform unimodular transformations in 
the space of the exponents to reduce the initial form system to a system of two polynomial in one 
variable. In [7], the coordinate transformations resulting from those unimodular transformations 
are called power transformations and they power up the field of "Power Geometry" . 

For example, investigating the direction (—1,-1): 

in ( _i -i)(/) = 54x 13 y 2 + 6x u y 

in ( _i _!)(#) = 72xV° + 8x 10 y 9 " 
We change coordinates using the unimodular matrix M = . 
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Definition 4.5 For a tropism (u,v) normalized so the greatest common divisor gcd(u,v) = 1, 
the unimodular matrix M 



M 



u v 
-I k 



?cd(u,v) = 1 = ku + lv = det(M) (11) 



defines the unimodular coordinate transformation x = X U Y 1 and y = X v Y k . 
Note that for a monomial x a y b , the coordinate transformation yields 

(X U Y — ') a (X V Y^ )k j^au+bvy — la+kb j^{(a,b),(u,v)) y — la+kb (12) 

so after the coordinate transformation, the monomials in the initial forms all have the same 
minimal degree in X. 

Continuing the example, we perform the change of coordinates: 

( m { _ h _ 1) (f)(x = X-\y = X- l Y- 1 ) = (54Y + 6)/(X 15 Y 2 ) 

< (13) 

I m { _ h _ 1) (g)(x = X- l ,y = X- l Y- 1 ) = (72Y + 8)/(X 19 Y w ) 

This change of coordinates reduces the initial form system to a system of two polynomials in one 
variable. For the example, Y = —1/9 represents the common root at infinity. Going back to the 
original coordinates: 

X = t ( x = X- x \ _ ( x = t- 1 

Y = -l/9 \y = X- 1 Y- 1 ) ^ \y = -9t-\ ^ 

As t goes to we have indeed a root going off to infinity. 

For every tentacle of the common factor we can associate a degree as follows. Considering 
again the common factor r from (JSj) , the amoeba for r has four tentacles, see Figure [7J reflected 
by its tropicalization 

Trop(r) = { (1, 0), (0, 1), (-1, -1), (0, -1) } . (15) 

In Table Q] we list the degrees associated to each tentacle of the common factor. We count the 
number of nonzero solutions of the initial forms, after proper unimodular coordinate transfor- 
mation. To make the correspondence with the usual degree, we consider only tropisms with 
nonnegative exponents and ignore zero factors. 

(u,v) in M (r) degree 

(1,0) 2xy + 9xy 2 1 

(0, 1) 2xy + x 2 y + 7x 3 y + x y 3 

(-1,-1) x 4 y + 9x 3 y 2 1 

(0,-1) 9xy 2 + 9x 3 y 2 2 

Table 1: Degrees associated to each vector in Trop(r). 



10 



4.3 A preprocessing algorithm and its cost 



That two polynomials with approximate coefficients have a common factor does not happen that 
often. Therefore, it is important to be able to decide quickly in case there is no common factor. 
The stages in a preprocessing algorithm are sketched in Figure 14.31 



tropicalization 

I 

1. compute tropisms O 



2. solve initial forms O 



3. compute 2nd term O 



no tropism 
=> no root at oo 

no root at oo 
no series 



no series 

no factor 



scries 



Figure 8: A staggered approach to computing a common factor. 

In Figure B~3l we distinguish three computational steps. We will address the cost of the first 
two steps in the following propositions. 

Proposition 4.6 Let f and g be two polynomials given by respectively n and m monomials. The 
cost of computing tropisms Trop(/) Pi Trop(g) is 0(nlog(n)) + 0(mlog(m)). 

Proof. It takes 0(nlog(n)) operations for computing a tropicalization Trop(/) because computing 
the convex hull of a set of n points amounts to sorting the points in the support. Likewise, 
computing Trop(g) takes 0(mlog(n)) operations. Merging sorted lists of normals to find the 
tropisms in Trop(/) H Trop(g) takes linear time in the length of the lists. □ 

This preprocessing step has the lowest complexity and as the algorithm operates only on the 
exponents the outcome is exact. The absence of tropisms is an exact certificate that there is no 
common factor, for any nonzero choice of the coefficients of the polynomials. 

In case we have tropisms, we have to compute roots at infinity. The cost of the second 
preprocessing stage is as follows. 

Proposition 4.7 Let f and g be two polynomials given by respectively n and m monomials. For 
every tropism t G Trop(/) n Trop(g) it takes at most 0((n + m) 3 ) operations to find a common 
solution in (C*) 2 to the initial form system defined by v. 
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Proof. For a tropism v, we solve the initial form system. In particular, an initial root z satisfies 



{ 



in v (/)(z)=0 
in v (s)(z) = 



z G (C*) 2 . 



(16) 



We perform a unimodular transformation so the tropism we consider is a unit vector, (1,0) or (0,1). 
This implies that the two equations in the initial form system are defined by two polynomials in 
one variable. To decide whether two polynomials in one variable admit a common solution we 
determine the rank of the Sylvester matrix. Using singular value decomposition, the cost of this 



Even as the cost estimates in the propositions are conservative, they give a good polynomial 
complexity. Actually, in the best case, the initial forms are supported on two points only and 
instead of a rank determination, we can just take primitive roots. The cost estimates of Propo- 
sition 14.71 cover the very worst situation where the Newton polygons are triangles and one of the 
edges contains all exponent vectors except for one. 

For numerical calculations, it is important to note that in this preprocessing stage, only the 
coefficients at the edges are involved. If the coefficients are badly scaled, then coefficients with 
monomials in the interior of the Newton polygons will not cause difficulties in this preprocessing 
stage. 

For the complexity in the proof of the second proposition we used the ubiquitous singular 
value decomposition but for practical purposes rank revealing algorithms [32J have a lower cost. 
The accurate location of the root of the initial form systems may look complication in case this 
root is multiple. However, because the initial form systems consists of univariate equations, the 
methods of [58] will give satisfactory answers. 

5 Germs 

Once we have a tropisms and an initial root at infinity, we start growing the Puiseux series for 
the common factor. 

5.1 how the amoeba grows from infinity 

We use the roots at infinity to grow the tentacles of the common factor. But first we must decide 
whether the roots at infinity are isolated or not. 

The canonical form of a fractional power series Puiseux series of a curve is defined next. 

Definition 5.1 Consider the curve defined by r(x,y) = 0. A fractional power series has the 
form 



The leading exponents (u,v) are the tropisms. The leading coefficients {zi,zz) satisfy the initial 



rank determination is cubic in the size of the matrix. 



□ 





(17) 
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form system 

™-(u !V )(f)(z 1 ,z 2 ) = 



m (u,v)(9)(zi,z 2 ) = ^ 

where the initial forms are taken from the equations / and g which define the common factor r. 

Continuing our example, we extend the solution at infinity, defined by the initial form system 
for the first tropism (1,0). Because the tropism is a standard basis vector, the Maple command 
sort({ f , g }, plex, ascending) will show that the leading terms of the polynomials / and 
g are indeed in (li0) (/) and m (1)0 )(#): 

{/ = 10xy 5 + 45xy 7 + 55xy 6 + x 2 { 30 other terms ) 
g = A5xy 7 + Wxy e + x 2 ( 34 other terms ) 
Let fi = f jx and g\ = g/x, then z = —2/9 is solution at infinity. 

y = -^t° + Ct(l + 0(t)), ceC*. (20) 

A nonzero value for C will give the third certificate of a common factor. Useful Maple commands 
to compute the power series are 

zt := x = t, y = -2/9 + C*t; 

flz := subs(zt,fl): glz := subs(zt,gl): 

cl := coef f (f lz,t , 1) ; c2 := coef f (glz.t , 1) ; 

The constraints on the coefficient C we obtain are 

1120 _ 1120 (~i _ n 
531441 59049° ~~ U . . 

320 _ 320 fi _ n ^ ' 

59049 531441 ° u 

Notice that the second coefficient C of the Puiseux series expansion again must satisfy an overde- 
termined system. Solving both equations for C gives C = —1/9. 



x = t 

y = -| - |*(1 + O(t)). 



(22) 



Substituting x = t,y = —2/9 — t/9 into f\ and g\ gives 0(t 2 ). The second term of the Puiseux 
series is the third and last certificate for a common factor. 

In general, the next term in the Puiseux series expansion might have a degree higher than 
one, or they might not exist a second term at all in case the solution at infinity is isolated. There 
is an explicit condition on the exponent of the second term in the Puiseux series expansion as in 
Proposition 15.21 

Proposition 5.2 Let f and g be two polynomials in two variables with (u, v) a tropism for which 
there is a root of the corresponding initial form system. After a unimodular coordination defined 
by the tropism (u, v) the Puiseux series takes the standard form 

n , n +w with d = Scd(u, v) (23) 
y — cq + c\t 
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and cq is a nonzero solution of the initial form system: 



(25) 



*f°) = n where ( ?W)tt x > V} = ^ (24) 
Assuming that 

f(x;y) =p{y) + x k P(x,y), k > 0, P(x,y) = p 00 + p w x + p 01 y -\ 

d{x,y) = q(y) +x l Q(x,y), I > 0, Q(x, y) = g o + qiox + goi2/ H 

so that we may write p and q as follows: 

V = 7o + 7iy ai + l2V a2 + ■■■ and q = S + 5 iy b i + 5 2 y b2 + ■■■ m) 
< a% < a 2 < ■ ■ ■ an < bi < b 2 < ■ ■ ■ 

then the condition on the exponent w such that there may exist a nonzero c\ is 

kd Id , , 

w = — = -. 27 
ai bi 



If*Mj:ix« } then there is no second term in the Puiseux series expansion and the root of the 
initial form system is isolated. 



Proof. We substitute the series (|23p into the assumed form of the system (J25J) and look for the 
lowest w so that the coefficient with the monomial of the lowest power of t vanishes. For the 
polynomials in the initial form system the substitution yields 

p(c + cit w ) = a lCl t aiW + h.o.t. and q(c + Cl t w ) = f3 lCl t blW + h.o.t. (28) 

and for the other terms x k P(x,y) and x l Q(x,y) we find 

t kd (a 2 + 0(t)),a 2 =R)o+Poic and t ld {[3 2 + 0{t)), fo = qoo + QoiCo- (29) 

Collecting the lowest degree terms in t we find the system 

aicit aiw + a 2 t kd = 



P lCl t b ^ w + f3 2 t ld = 0. ( 3 °) 

For ci to exists as a nonzero solution of a linear system annihilating the lowest degree terms in t, 
the exponents of t in the system (|30j) must be equal in both equations: 

a\W = kd and b\w = Id. (31) 

Eliminating w leads to the condition (|27p . □ 



The assumption in <\25\\ of Proposition 15.21 makes abstraction of the monomial shifts that are 
usually required when looking for nonzero solutions. So the assumption (|25p will not hold in 
general and monomial shifts are needed. But this will still lead to an explicit condition on the 
exponents purely derived from the exponents of the polynomials / and g for a second term in the 
Puiseux series to exist, leading to a linear system of two equations in the one coefficient c\. 

In case the condition of Proposition 15.21 is satisfied and the linear system admits a nonzero 
solution for ci, then the exponent w and coefficient c\ constitute respectively an exact and an 
approximate certificate for the existence of a common factor for the two polynomials / and g. 
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5.2 regions of convergence of Puiseux series 

Neighborhoods where Puiseux series converge are called germs. Following }14j , we get the defini- 
tion for a germ: 

Definition 5.3 Given a power series expansion at a point, a germ is an equivalence class of open 
neighborhoods of the point where the series converge. 

Via substitution as in the proof of Proposition 15.21 we may compute more terms in the Puiseux 
series. In each step, the exponent in t will increase by at least one. In each step, also more and 
more approximate coefficients of the given polynomials will be used. The accuracy of the given 
input coefficients will impose a natural bound on the order for which it still makes sense to extend 
the series. 

The unimodular coordinate transformations play a very important role also in the accurate 
evaluation of polynomials [15J . As the size of arguments of the polynomial functions grows, and as 
the direction of the growth points along the direction of a tentacle of the amoeba, monomials on 
the faces perpendicular to that direction become dominant. A weighted projective transformation 
as in [53j will rescale the problem of evaluating a high degree polynomial with approximate 
coefficients near a root. 

6 Implementation Aspects 

For efficient implementation of the algorithm, the data structures used to represent the polyno- 
mials consist of a list of exponent vectors and a coefficient table. More precisely, to represent a 
polynomial / denoted as 



we use a list to represent the support A and a lookup table C[A] for the coefficients. The indices 
of the lookup table Ca are the exponent vectors a € A. In Maple's index notation: C[a] = c a . 

Separating the support from the coefficient allows an efficient execution of change of monomial 
orders. If n = #A, then monomial orders on / are stored via permutations of the first n natural 
numbers. The separation also gives an efficient way to change coordinates, i.e.: we apply the 
unimodular coordinate transformation only on A. For a unimodular matrix M: 



Abusing notation, for zEC: Mz denotes the value for Y after applying the coordinate trans- 
formation as in Definition 14.51 

The input polynomials / and g with respective supports Af and A g are then represented by 
two tuples: (Af,C[Af]) and (A g ,C[A g ]). The preprocessing algorithms consists of two stage. In 
the first stage, Algorithm 16.11 computes the tropisms and the roots of the corresponding initial 
form systems. If the sets of roots are not empty, the exponent and coefficients of the second term 
in the Puiseux expansions are computed by Algorithm 16.21 in the second stage. We define the 
specifications of the algorithms below. 




(32) 



aeA 



MA = { Ma j a e A }. 



(33) 
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Algorithm 6.1 Tropisms and Initial Roots 

Input : (Af,C[Af)) and (A g ,C[A g ]). 
Output : T = { (u,v) (£Z 2 \ (0, 0) | («, v) is tropism }, 

R[T] = { { z € C* | in (U)U) (/)(M*) = 0,in M (j)(M 2 ) = } | (u, «) G T }. 

Every tropism in T defines a set of roots (possibly empty) of the corresponding initial form system, 
after application of the unimodular coordinate transformation M. The cost of Algorithm 16,11 is 
estimated by Proposition 14.61 and Proposition 14.71 

Algorithm 6.2 Second Term of Puiseux Expansion 

Input : (Af, C[A f ]), (A g , C[A g \), T, and R[T}. 
Output : W[i2[T]] = { (c, w) G C* x N+ | z G Z G R[T] }. 

The elements of the set W[i?[T]] define the second term of the Puiseux series expansion. In 
particular, for every (c,w) G W[i2[T]]: 

f X = id 

where (X, Y) are the new coordinates after applying the transformation of Definition 14.51 Con- 
ditions on the existence of the exponent w are given in Proposition 15.21 

7 Conclusions and Extensions 

Like Maple, tropical algebraic geometry is language. Sentences like tropisms give the germs to 
grow the tentacles of the common amoeba express efficient preprocessing stages to detect and 
compute common factors of two polynomials with approximate coefficients. Seeing the problem 
as a system of two polynomial equations in two variables, the algorithm is a polyhedral method 
to find algebraic curves. 

Among the extensions we consider for future developments are algorithms to handle singular- 
ities numerically and polyhedral methods for space curves. 
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